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WAVE  PROPAGATION  IN  LAMINATED  COMPOSITE  PLATES 


Introduction 


In  the  past  dynamic  behavior  of  infinite  periodically  laminated  medium  has  been  stu* 
died  extensively.  A  review  of  the  literature  on  exact  and  approximate  analysis  of  this 
problem  can  be  found  in  [1.2].  In  [2]  a  stiffness  method  was  presented  for  studying 
harmonic  wave  propagation  in  a  periodically  laminated  infinite  medium.  In  this  method 
each  lamina  is  divided  into  several  sublayers  and  the  displacement  distribution  through  the 
thickness  of  each  sublayer  is  approximated  by  polynomial  interpolation  functions.  These 
functions  involve  number  of  discrete  generalized  coordinates,  which  are  the  displacements 
and  tractions  at  the  interfaces  between  the  adjoining  sublayers.  This  ensures  continuity  of 
these  quantities  at  these  interfaces.  By  applying  Hamilton’s  principle  the  dispersion  equa¬ 
tion  is  obtained  as  a  standard  algebraic  eigenvalue  problem.  Eignevalues  and  eigenvectors 
of  this  equation  yield  the  frequencies  of  propagating  and  evanescent  modes  and  the  associ¬ 
ated  displacements  and  tractions  at  the  intefaces.  Results  obtained  by  this  method  for  an 
infinite  medium  with  isotropic  laminae  were  shown  to  agree  well  with  exact  solutions  [1]. 


In  this  paper  we  use  a  similiar  technique  to  study  dispersion  of  waves  in  a  lami¬ 
nated  plate.  Although  the  method  is  easily  applied  to  arbitrarily  anisotropic  laminae,  we 
have  illustrated  the  method  for  the  case  when  each  lamina  is  transversely  isotropic  with 
the  symmetry  axis  lying  parallel  to  the  lamina.  This  study  is  motivated  by  the  desire  to 
model  wave  propagation  in  a  continuous  fiber-reinforced  laminated  plate.  If  it  is  assumed 

that  the  wave  lengths  are  much  longer  than  the  fiber  diameters  and  spacing,  then  each 

lamina  can  be  modeled  as  a  homogeneous  transversely  isotropic  medium  with  the  symmetry 
axis  paralled  to  the  fibers.  The  overall  effective  elastic  properties  of  such  a  medium  can 
be  calculated  from  the  fiber  and  matrix  properties  by  using  an  effective  modulus  theory 
[3,4].  Such  an  assumption  has  been  made  in  this  paper.  Thus  this  study  examines  the 
dispersion  characteristics  of  guided  waves  in  a  layered  anisotropic  plate. 

Although  wave  propagation  in  a  homogeneous  isotropic  plate  has  been  thoroughly 
studied  using  exact  or  approximate  means  by  Mindlin  [5]  and  many  others,  not  many 

exact  solutions  for  layered  plates  have  been  reported.  Most  of  the  works  on  layered 


plates  have  used  approximate  plate  theories  or  effective  modulus  theory  for  the  entire  plate 
(Dong  and  Nelson  [6]  for  discussion)  because  of  the  extreme  complexities  involved  in  the 
equations  obtained  when  the  exact  three  dimensional  linear  theory  of  elasticity  is  used. 
For  a  two  layer  or  three  layer  (sandwich)  isotropic  plate  the  governing  dispersion  equation 
becomes  tractable  and  has  been  discussed  in  [7-9]. 

Wave  propagation  in  a  monoclinic  crystal  plate  has  been  fully*- investigated  in  [10-12]. 
However,  to  our  knowledge,  no  exact  solutions  for  laminated  anisotropic  plate  are  avail¬ 
able.  In  [6.13.15]  the  authors  present  a  stiffness  method  to  study  this  problem.  They 

use  interpolation  functions  for  the  displacements  that  involve  as  generalized  coordinates  the 

displacements  at  the  interfaces  between  layers  and  at  the  middle  of  the  layers.  Although 

>  f 

the  technique  yields  generally  good  results  for  the  frequency  spectrum,  it  does  not  lead 
directly  to  the  determination  of  the  tractions  at  the  interfaces. 

For  studying  scattering  by  cracks  or  inhomogeneit.es  using  a  hybrid  method  in  which 
the  localized  region  near  the  scatterers  is  represented  by  finite  elements  and  the  exterior 
solution  by  superposition  of  modes  (propagating  and  non-propagating)  it  is  necessary  not 
only  to  get  accurate  approximations  to  displacements,  but  also  to  the  tractions  at  the 

boundary  separating  the  inner  and  outer  regions  (for  applications  of  this  hybrid  technique 

see  [15-17]).  It  is  with  this  application  in  mind  that  we  present  in  this  paper  a  stiffness 

method  that  incorporates  as  generalized  coordinates  displacements  and  tractions  at  the  inter¬ 
faces  between  the  layers  (and  sublayers).  As  mentioned  before  the  results  obtained  by 
this  method  for  a  periodically  laminated  infinite  medium  agree  well  with  known  exact  solu¬ 
tions.  Here  we  show  that  for  an  isotropic  and  transversely  isotropic  plate  the  method 
yields  results  in  excellent  agreement  with  the  exact  solutions.  It  is  noted  that  the  present 
method  using  the  same  number  of  sublayers  yields  more  accurate  results  at  high  frequencies 
than  obtained  by  the  technique  of  Dong  et  al  [6,13,14].  Intuitively  it  is  clear  that  as  the 
number  of  sublayers  tend  to  infinity  the  results  based  on  the  present  formulation,  where 
the  continuity  of  displacement  and  traction  is  maintained  at  the  interfaces,  should  approach 
the  exact  continuum  solution. 
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Governing  Equations 


We  consider  time  harmonic  waves  in  a  plate  composed  of  n  laminae.  For  simplicity 
in  analysis  it  will  be  assumed  that  each  lamina  is  transversely  isotropic  with  the  symmetry 
axis  aligned  with  either  the  x  or  the  y  axis  (Fig.  1).  This  assumption  is  not  necessary 
for  the  development  of  the  equations  but.  it  is  made  here  to  keep  the  algebra  as  simple  as 

-V 

possible  for  the  anisotropic  problem  at  hand.  Under  this  assumption  the  wave  propaga¬ 
tion  problem  reduces  to  two  uncoupled  ones:  plane  strain  in  which  the  displacement 
components  are  ux.  o,  uz,  and  SH  or  antiplane  strain  when  the  only  nonzero  displace¬ 
ment  component  is  uy.  So  in  the  following  we  shall  treat  these  two  separately. 


Plane  Strain 


Consider  the  ith  lamina  bounded  by  z  -  zj.j  and  z  -  zj+j.  The  stress  strain 
relation  in  this  lamina  will  be  given  by 


c(i)  c(i)  0 
11  13 

c®  c(i)  0 
13  33 
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where  ajj  and  ejj  are  the  stress  and  strain  components,  respectively,  and  we  have  written 
7Xz  "  2«xz.  Note  that  if  the  y  axis  is  the  axis  of  symmetry,  then 


(i)  -  c(i)  and  c(i)  -  *fc(i)  -  c(i)\ 
11  33  55  2\  11  13  j 
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Then  the  problem  is  equivalent  to  that  of  an  isotropic  one.  In  order  to  get  good  numer¬ 
ical  results  each  lamina  will  be  divided  into  several  sublayers,  m.  say.  Within  the  ith 
sublayer  we  will  choose  a  local  coordinate  with  the  origin  at  the  mid-plane  and  Xj.  yj. 
zj.  parallel  to  the  global  x.y.z  axes, 

respectively.  Let  2hj  be  the  thickness  of  this  sublayer.  Denoting  u(i>  to  be  the  displace¬ 
ment  at  a  point  in  the  jth  lamina  we  write  ■> 
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where  fn  (n-1,  ....  4)  are  cubic  polynomials  in  the  local  coordinate  Zj  given  by 


f,  -  i  (2-3nj  +  n’j).  ft  “  \  <2+3i|j  -  rj’j) 


f,  -  -j-  (1-Oj  -  V’j  *  v’j).  U  -  -j-  (-1  -»»j  ♦  n*j  +  n’j) 


Here  -  zj/hj  and  uj.  wj,  Xj,  ej  are  the  values  of  ux,  uz.  <rX2.  and  ozz  at  the  jth 

node.  These  nodal  values  of  the  displacement  and  traction  components  are  functions  of 
xj  (-x)  and  t.  In  this  paper  it  will  be  assumed  that  the  time  dependence  is  of  the  form 

e-iuft  gi>  being  the  circular  frequency.  The  factor  e"*®1  will  be  dropped  in  the  sequel. 
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The  equations  governing  the  nodal  generalized  coordinates  {uj.  Xj .  wj.  ?.}  will  be 
obtained  using  Hamilton's  principle.  For  this  purpose  we  calculate  the  lagrangian,  L.  per 
unit  length  in  the  y-direction  of  the  plate  as 


t  >■»  -th'  •  0«1>)  Ur* 


i-1  j-1 


•  L  I  «■». 


i-l  j-1 


The  overbar  denotes  complex  conjugate. 


Using  (2)  in  the  strain-displacement  relation,  and  in  turn  in  (4),  we  get  the  expres¬ 


sion  for  L*‘)j  as 
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where  {q}  is  defined  as 


S'  "^J-4  x 

(q)  -  £uj-i.  Xjt{.  wj.j,  0j.x.  Uj.  xj.  wj.  ajj  T 


and  the  primes  denote  derivatives  with  respect  to  x.  The  matrices  [C,].  [C,].  etc.  are 
defined  in  the  Appendix  1. 


The  Lagrangian  for  the  complete  plate  is  obtained  by  summing  over  all  the  layers 
and  its  variation  leads  to  the  governing  equation  for  the  plate  as 


[Q"]  -  [CJ  [Q’]+  [M]  [Q]) 


-  (fEJ  CQiv3  ♦  [EJ  [Q***]  ♦  [E,}  [Q'*3  ♦ 


[EJ  fQ  ]  +  IE,]  [Q])  -  0 


Note  that  the  matrices  [C,],  [M],  [E,],  [E,],  and  [E,]  are  symmetric,  whereas,  the  others 
are  skew-symmetric.  These  assembled  matrices  are  defined  in  Appendix  I.  It  is  seen  that 
the  generalized  coordinate  vector  {Q}  satisfies  a  fourth  order  homogeneous  ordinary  differ¬ 
ential  equation  in  x. 

A  solution  to  equation  (7)  can  be  assumed  in  the  form 


tQ]  -  [Q.]e*x 


where  [QJ  represents  the  amplitude  vc-tor.  Substituting  (8)  in  (7)  we  get  a  set  of  linear 
homogeneous  equations  to  solve  for  [QJ: 


k4[K,]  -  ik’[KJ  -  k*[KJ  +  ikCK.]  ♦  [K,]  {Q.}  -  0 


(9) 


where 

[KJ  -  [EJ 

[KJ  -  [EJ 

%  ' 

[KJ  -  [EJ  +  u’[CJ 

[KJ  -  [EJ  +  «*[CJ 

[KJ  -  [E,j  -  «*[M] 

For  nontrivial  solution  the  determinant  of  the  coefficient  matrix  must  be  zero.  This 
equation  is  the  dispersion  equation  to  solve  for  the  eigenvalues  k  for  given  w.  Alterna¬ 
tively,  Eq.  (9)  can  also  be  written  as  ([Kp]  -  w’[MpJ)  [QJ  -  0 
where 


[Kp]  -  k4tKJ  -  ik*[KJ  -  k’[EJ  ♦  ik[EJ  +  [EJ,  [Mp]  -  [M]  -  ik[CJ  +  k’[CJ 


This  leads  to  the  eigenvalue  problem  for  solving  u1  for  given  k. 


Antiplane  Strain 


The  derivation  given  above  is  for  the  case  of  plane  strain.  The  antiplane  strain 
case  can  be  considered  in  a  similiar  manner  and  that  results  in  a  much  simpler  fre¬ 
quency/wave  number  equation.  In  this  case  we  assume 


V:  i  f.  +  V:  f,  +  ~ —  r:i  fj  +  — - —  r:  f4 
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The  corresponding  eigenvalue  problem  for  the  wave  number  is 

<-k*[KJ  +  CK.])  {Q.}  -  0 

and  that  for  the  frequency  is 

(CK,]  -  cj*[Ms])  {Q.}  -  0 

where 


CK,]  -  kJ[KJ  ♦  [Kt] 

CK,]  -  [Ms]  -  [K,] 


(12) 


and  the  corresponding  matrices  for  the  sublayer  from  which  these  are  assembled  are  def- 


ined  in  Appendix  I.  The  vector  {q}  is  defined  as  Jvj.j.  rj.j.  Vj,  r j  | 

As  illustrative  examples  we  show  equations  (9)  and  (10)  for  an  isotropic  plate,  a 
fiber-reinforced  ((f)  homogeneous  plate,  a  three-layered  (0*/90*/0*)  plate  and  a  four-layered 
(0*/90,/0*/90*)  plate.  These  results  are  discussed  in  the  next  section-  We  also  present 
some  numerical  results  for  the  SH  case. 

Numerical  Results  and  Discussion 

In  orde'  to  validate  the  method  we  first  considered  an  isotrt>pic  plate  with  Poisson's 
ratio,  r  -  0.31.  A  full  discussion  of  the  frequency  spectrum  for  this  case  was  given  by 
Mindlin  [5].  Both  the  eigenvalue  problems  (9)  and  (10)  were  solved.  The  comparison  of 
the  results  obtained  for  real  and  imaginary  wave  numbers  with  those  of  Ref.  [5]  are 

shown  in  Fig.  2.  In  this  figure  8  -  - .  y  ■  — .  where  H  is  thickness  of  the 

wzr, 

plate.  For  this  computation  the  plate  was  divided  into  25  sublayers.  It  is  seen  that  the 
results  agree  with  the  exact  solution  extremely  well  at  low  and  moderate  frequencies.  The 
agreement  at  high  frequencies  is  also  quite  reasonable.  We  also  used  the  interpolation 
functions  adopted  by  Dong  and  Huang  [15],  and  obtained  the  frequency  spectrum  keeping 
the  same  number  of  sublayers.  Comparison  of  the  results  with  the  exact  solution  is 

shown  in  Fig.  3.  The  loss  of  accuracy  at  high  frequencies  is  evident.  Of  course,  the 

accuracy  will  improve  if  more '  sublayers  are  taken  [6], 

We  then  considered  a  fiber-reinforced  plate  when  the  fibrs  are  aligned  along  the  x- 
axis  (0*) .  The  properties  of  the  plate  are  given  in  Table  1.  As  mentioned  before,  for 

propagation  in  the  0*  (90*)  direction  the  problem  is  tractable  anaiytically.  The  frequency 

equation  for  this  case  is  given  in  Appendix  II.  Figure  4  shows  the  comparison  of  the 
frequency  spectrum  obtained  by  the  present  method  with  the  exact  solution  for  propagation 


in  the  x-direction.  it  is  seen  that  the  comparison  is  excellent.  In  this  figure 


kH 

-  —  and  7  *  — .  It  is  seen  that  as  k  «  the  slopes  of  the  disper- 
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sion  curves  for  the  first  symmetric  and  antisymmetric  modes  tend  to  the  ratio 

VR/JC55/P.  where  Vj^  is  the  Rayleigh  wave  velocity  in  the  x-direction.  For  the  dis¬ 
cretization  method  we  used  15  sublayers. 

The  next  problem  considered  was  a  sandwich  plate  with  0°/90°/0°  configuration. 
Figure  5  shows  the  dispersion  curves.  These  curves  differ  considerably  from  those  shown 
in  Fig.  4.  It  is  seen  that  the  cut-off  frequencies  are  lowered  as  jvell  as  the  slopes  of  the 
curves.  Branches  of  imaginary  wave  numbers  are  also  quite  different.  In  these  computa¬ 
tions  we  used  15  sublayers  also. 

To  see  the  comparison  between  our  results  and  those  obtained  by  the  author  of 

Ref.  [14]  we  considered  a  4-layer  (0*/90*/0“/90“>  plate  with  the  properties*  given  in  [14]. 

Solid  lines  in  Fig.  <5-  are  Dong  and  Huang’s  results  using  20  sublayers.  Open  circles  and 

dashed  lines  are  ours  using  the  same  number  of  sublayers.  The  agreement  deteriorates  at 
high  frequencies.  For  the  same  composite  plate  we  show  in  Fig.  7  the  complex  frequency 
branches.  Solid  dots  are  obtained  by  the  method  of  [14], 


Finally,  figures  8  and  9  show  the  SH  wave  spectra  for  homogeneous  (0°)  and  three- 
layer  (0*/90“/0*)  plates.  Although  the  curves  look  very  similiar,  it  is  interesting  to  note 
that  the  cut-off  frequencies  are  increased  in  the  latter  case.  This  is  contrary  to  what  was 
found  in  the  plane  strain  case. 


*(Cu)o*  ’  21.289  X  10*  psi;  (C„)(j.  -  0.592  X  10*  psi;  (C„)o»  -  2.3186  X  10*  psi;  (C„)q» 


■  0.850  X  10*  psi;  p  *  1  -®-r 

cm 
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Conclusion 


A  stiffness  method  using  through- thickness  interpolation  functions  for  the  displace- 
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ments  that  maintain  continuity  of  displacements  and  tractions  at  the  interfaces  between  the 

layers  has  been  used  to  study  guided  waves  in  a  laminated  composite  plate.  It  is  shown 

that  the  results  agree  with  exact  known  solutions  for  homogeneous  plates  and  also  with 

their  numerical  solutions.  The  advantage  of  the  method  is  that  it  is  more  accurate  at  high 

frequencies  with  smaller  number  of  sublayers  than  others  in  which  only  displacement  con- 

»  1 

tinuity  is  enforced.  Also,  tractions  and  displacements  are  obtained  as  eigenvectors,  so  that 
the  former  need  not  be  calculated  by  differentation  afterwards. 

Numerical  results  presented  show  significant  differences  between  the  dispersion  curves 
for  the  homogeneous  fiber-reinforced  plate  and  the  cross-ply  laminated  plate. 

Although  the  case  of  wave  propagation  along  a  principal  direction  has  been  studied, 
the  method  is  easily  applicable  to  off-axis  propagation.  For  purposes  of  illustration  only 
a  small  number  of  laminae  has  been  considered.  But,  because  we  divide  each  lamina  into 
several  sublayers,  each  of  these  sublayers  can  have  different  properties.  So  any  number  of 
laminae  can  be  considered. 
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Table  I.  Properties  of  0*  and  90*  laminae.  All  the  stiffnesses  are  in  units  of  10"  N/mJ. 


pfg/cm1)  C„  Cjj  C  |2  C»4  CM 

0*  lamina  1.2  1.6073  0.1392  0.0644  0.0350  0.0707 

90*  lamina  1.2  0.1392  1.6073  0.0644  0.0707  0.0350 
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Appendix  I 


We  first  define  the  following  matrices: 
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[CJ,  etc.  can  be  written  as 
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For  antiplane  motion  we  define: 
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Appendix  II 


For  the  propagation  of  waves  along  the  plate  of  transversely  isotropic  material 
considered  here  the  dispersion  equation  separates  into  two.  one  for  the  symmetric  motion 
and  the  other  for  the  antisymmetric  case.  In  the  former  case  we  get 

tan(s,H/2)  s,E,T, 

tan(s,H/2)  “  s,E,T,  111 ’ 11 

whereas,  for  the  latter  we  have 

tan(s,H/2)  _  s^T, 

tan(s,H/2)  ’  stE,T,  '  2) 


In  writing  (II. 1)  and  (11.2)  we  have  defined  the  following  quantities, 
a  *•  C|,/Cu<  P  m  Cjj/Cu 
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